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Abstract. We study an example of instability in presence of a multiplicative noise, namely the spontaneous 
generation of a magnetic field in a turbulent medium. This so-called turbulent dynamo problem remains 
challenging, experimentally and theoretically. In this field, the prevailing theory is the Mean-Field Dynamo 
[T] where the dynamo effect is monitored by the mean magnetic field (other possible choices would be the 
energy, fiux,...). In recent years, it has been shown on stochastic oscillators that this type of approach 
could be misleading. In this paper, we develop a stochastic description of the turbulent dynamo effect 
which permits us to define unambiguously a threshold for the dynamo effect, namely by globally analyzing 
the probability density function of the magnetic field instead of a given moment. 

PACS. 02.50.-r Probability theory, stochastic processes and statistics - 47.27.Gs Isotropic turbulence; 
homogeneous turbulence - 47.27.Jv High-Reynolds-number turbulence 

Classical stability analysis arc usually performed in idenced on a DufRng oscillator , where the solution 

systems were the control parameter is a non-fluctuating x(t) ~ is stable for values of the control parameter above 

quantity, e.g. for laminar flows in hydrodynamics. When the deterministic threshold. This stabilization is generic 

the instability occurs in a random system (e.g. a turbu- for weak intensities of the noise. For stronger intensity, 

lent medium), resulting fluctuation of the control param- however, noise induced transition may also arise in this 

eter, or multiplicative noise, may generate several surpris- system [3]. In the case of a parametric instability, it has 

ing effects that have been studied in a variety of systems, been showed experimentally 4J that the instability is sen- 

The possibility of stabilization by noise has flrst been ev- sitive to the bifurcation nature. In the supercritical case. 
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oscillatory bursts (corresponding to a signal with a vanish- 
ing mean but a most probable value equal to zero) appear 
first when the control parameter is increased and are then 
replaced by a state where the most probable value is no 
more equal to zero. On the contrary, in the subcritical 
case, there is coexistence between these two states. This 
illustrates a central difficulty of instability in presence of 
multiplicative noise associated with an ambiguity regard- 
ing the threshold value, which depends on the definition 
of the order parameter |5] . 

The observations and techniques developed in these 
simple systems may be used to shed new light on some 
recent issues associated with the dynamo effect, the pro- 
cess of magnetic field generation through the movement 
of an electrically conducting medium. In this case, the in- 
stability results from a competition between amplification 
of a seed magnetic field via stretching and folding, and 
magnetic field damping through diffusion. In a laminar 
fluid, it is controlled by a dimensionless number, the mag- 
netic Reynolds number (Rm), which must exceed some 
critical value Rrric for the instability to operate. In a tur- 
bulent medium, velocity fluctuations induce fluctuation 
of the control parameter, making the turbulent dynamo 
problem similar to an instability in the presence of multi- 
plicative noise. In that respect, recent numerical findings 
such as observed in [6j may find a natural explanation. In 
their work, the authors observed short intermittent bursts 
of magnetic activity separated by relatively long periods, 
increasing towards the bifurcation threshold. This feature 
could be explained in terms of a supercritical instability in 
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presence of multiplicative noise since in this case, the bi- 
furcated state is generally composed of oscillatory bursts. 
More generally, the multiplicative noise paradigm could 
turn useful to interprete the outcome of recent experi- 
ments involving liquid metals. Among the various operat- 
ing experiments, a clear distinction appears between set 
up with constrained or unconstrained geometry. In the for- 
mer case [7 , the fluctuation level is very weak. The veloc- 
ity field is then very close to its laminar (mean) value. In 
these experiments, dynamos have been observed, at crit- 
ical magnetic Reynolds number comparable to the theo- 
retical value. In contrast, unconstrained experiments [5] 
are characterized by a large fluctuation level (as high as 
50 per cent). A surrogate laminar Rrric can then be com- 
puted, using the mean velocity field as an input [9 but it is 
not clear whether it will correspond to the actual dynamo 
threshold, owing to the influence of turbulent fluctuations. 

In this paper, we investigate this issue by techniques 
developed to study the Duffing oscillator and using a stochas- 
tic description of small-scale turbulent motions. This sub- 
ject has been pioneered by Kazantsev [TD], Parker |TT] 
and Kraichnan 12J, and further developed by the Russian 
school [13]. It has recently been the subject of a renewed 
interest, in the framework of anomalous scaling and in- 
termittency [14j . or computation of turbulent transport 
coefficients and probability density functions (PDF) [15] . 

The dynamic of the magnetic field B in an infinite 
conducting medium of diffusivity rj and velocity V, is gov- 
erned by the induction equation: 

dtB, - -VkdkB, + BkdkV, + ridkduB, , (1) 
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with control parameter built using typical velocity and 
scale as Rm — LV/rj. We decompose the velocity field 
into a mean part Vi and a fluctuating part Vi. In most 
laboratory experiments, the mean part is provided by the 
forcing. As such, it is generally composed of large scales, 
while the fluctuating part collects all short time scale, 
small-scale movements. In this regard, it is natural to con- 
sider the fluctuating part of the velocity as a noise, to be 
prescribed or computed in a physically plausible manner. 
The simplest, most widely used shape is the Gaussian, 
delta-correlated fluctuations, the so-called "Kraichnan's 
ensemble" : 

{v,i^,t)v,i^',t')) = 2G,,(x,x')<5(i-t') (2) 

Equation ^ then takes the shape of a stochastic par- 
tial differential equation for B. In that respect, we note 
that the induction equation is linear and does not include 
explicit back reaction term allowing saturation of any po- 
tential growth in the dynamo regime. This back reaction 
is provided through the velocity which is subject to the 
Lorentz- Force, a quadratic form of _B. It is usually ignored 
in the so-called kinematic regime. However, for reasons 
which will become clearer later, we prefer to work with 
a modified induction equation, so as to model this non- 
linear back reaction. 

Indeed, the study of the induction equation alone could 
be misleading when looking at the threshold of the dy- 
namo effect (it leads to a threshold value dependent on 
the considered moment). This is symptomatic of the fact 
that the dynamo problem is nonlinear through the Navier- 
Stokes equation. A practical way to include the effect of 
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the Lorentz force at the onset of the nonlinear regime is 
to add a saturating term in the induction equation. Sym- 
metry considerations then favor a term like —cB^Bi. In 
some sense, this modification is akin to an amplitude equa- 
tion, and the cubic shape for the non-linear term could be 
viewed as the only one allowed by the symmetries. Such 
a procedure has been validated by [IS] in the case of the 
saturation of a Ponomarenko dynamo. Such a cubic form 
has also been evidenced by Boldyrev [15j by assuming the 
equality of viscous and dynamical stresses in the Navier- 
Stokes equation at the onset of backreaction. In the sequel, 
we show that the precise form of the nonlinear term does 
not affect the threshold value, which only depends on the 
behaviour for \B\ 0. This is similar to the case of the 
Duffing oscillator where the threshold is given by the Lya- 
punov exponent of the linearized oscillator ^ . 

A further difficulty is associated with the presence of 
the diffusive terms. Their physical influence is to cut-out 
magnetic field Fourier components with wave numbers ini- 
tially oriented in the contracting directions (stretching di- 
rections for wave numbers). Through the divergence free 
condition, this favors magnetic field components point- 
ing towards contracting directions, thus counteracting the 
initial effect of magnetic growth along stretching direc- 
tion. The balance between these two effects yields the dy- 
namo threshold ([H])- The direct consideration of diffu- 
sive terms in the stochastic formulation requires functional 
derivative and integration, hindering simple analytical de- 
scription. We propose to model them partially through an 
additional "molecular" homogeneous noise S,i(t), super- 
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posed to and decorrelated from the velocity fluctuation the mean-field, while the second one is akin to a turbulent 

|17l , with correlation function < >— '2.ri5ij5{t — diffusivity. A few remarks are in order at this point: i) our 

t'). In the sequel, we show that this choice provides some mean field equation has been derived without assumption 

sort of saturation for the moments of various order, sim- of scale separation, ii) The tensor ^ does not appear at 

ilar to the role of a viscosity. Damping of magnetic field this level. In the sequel, it will be shown to govern the 

fluctuations, however, is not properly taken account by stochastic dynamo transition, 
this model. 

Using standard techniques [l9lfT5]. one can then derive For this, we need to identify the threshold as a func- 

the evolution equation for P(B,x, t), the probabifity of tion of the noise properties. Here, we follow an idea by 

having the fleld B at point x and time t (we assume an Mallick and Marcq [.3_ , and focus on the properties of the 

homogeneous turbulence for simplicity): stationary PDF of the system. Indeed, below the transi- 
tion, the only stable state is S = and the PDF should be 



dtP = -VAP - {dkV,)dB, [BkP] + dk [hidiP] (3) 
+ c9b. [B''B,P] + 2dB, [BkaukdiP] 



a Dirac delta function. Above the transition, other equi- 
librium states are possible, with non zero magnetic field. 
fJ-tjkidBi[BjdBk{BiP)] , However, in the general case, it is not possible to find 

analytical solution for the equation We thus resort 
to the following mean field argument. Changing variable 
Bi = Bci where e is a unit vector (and can be character- 
and iJijki — {djVidiVk) ■ ized by d — 1 angular variables), we can get an equation 

for P{B,e^,x) = JP{B,) where J = B'^-^ is the Jaco- 
bian of the transformation. We now assume that there is 
an uncoupling for P as P{B^ ei) = P{B)G{ei, x), and per- 
form an average over the angular variables, to find a closed 
equation for P{B). In some sense, this can be regarded as 
a kind of angle/action variable separation, with average 
dt{B^) = -Vkdk{B') + {dkVt){Bk) - 2aktidk{Bi) q^^^. ^j^g f^st variables. The final equation for P{B) be- 

+ PkidkdiiB,) - c{B^B,). (5) comes: 

This equation resembles the classical Mean Field Equation 
of dynamo theory, with generalized (anisotropic) "a" and 

dP d d d d 

"/3". The first effect leads to a large scale instability for ^ = "'QB^^dB^^^^^ ^ ^95^^^^ ^ ""dB^^^^^ 



with the following turbulent tensors: 

f^ki = {vkVi) +r]Sku ciijk = {vidkVj) (4) 



The physical meaning of these tensors can be found 
by analogy with the "Mean-Field Dynamo theory" [T1[5D] . 
Indeed, consider the equation for the evolution of the mean 
field, obtained by multiplication of equation ^ by and 
integration: 
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where the coefBcients are given by averages over the posi- In some sense, this bifurcation (I) is obtained using the 

tion and the angular variables (m)^ — J • G{e^ x) dx.de: mean field as control parameter. Another bifurcation thrcsh- 

, , ,_^ old can be defined using the most probable field as control 

a = {fJ-ijkieiejekei}^ (7) 

parameter. Indeed, an elementary calculation shows that 

b = {dkVieiek) 4, + {fiijkii^ikejei + AkjCiei))^ 

the condition for a maximum in the PDF is 6 > a. There- 

where we used Z\y = dmirij) = Sij — CiCj an "angu- 

fore, the bifurcation threshold (II) with the most probable 

lar Dirac tensor" . One can notice that these coefficients 

field as control parameter is defined hy b = a. This differ- 
only explicitly involve the tensor fi. Nevertheless, one must 

ence may have some relevance when analyzing real data 

keep in mind that the tensor a and /3 enter these expres- 

from experiment. To illustrate this, we performed simula- 

sions by mean of the angular distribution G'(ei,x), whose 

tions of the 1-D version of our non-linear stochastic sys- 

expression involves this two tensors in the general case. 

tem. It may be checked that the stationary PDF in this 

An obvious stationary solution of ([5]) is a Dirac func- 

case is exactly given by eq. Therefore, the time series 

tion, representing a solution with vanishing magnetic field. 

and associated PDF are good illustration of the output of 

Another stationary solution can be found by setting = 

our 3D model, and the meaning of the two bifurcations 

in ([S]), with solution: 

discussed above. The time series and PDF for three dif- 

1 r c 1 

-P(^) = ^ ~ fcrcnt values of the control parameter are shown on figure 

where Z is a normalization constant. This solution can E 
represent a meaningful probability density function only 

an intermittent behavior for the magnetic field rcminis- 



The simulations show that the bifurcation (I) leads to 

if it can be normalized. This remark provides us with a 



bifurcation threshold: there is dynamo whenever ([81) is in- 

cent of the characteristic behavior of instability in pres- 

tegrable, i.e., when solution other than vanishing magnetic 

ence of multiplicative noise: typically, the magnetic energy 

field are possible. 

exhibits bursts separated by long quiescent (zero magnetic 

Condition of integrability at infinity of ^ requires a . . 

energy) period. Equation ([8]) cannot rigorously capture 

be positive. This illustrates the importance of the non- 

these intermittent states. However, two facts are very sug- 

linear term which is essential to ensure vanishing of the 

gestive of such a type of bifurcation in our solution: (a) the 

probability density at infinity. Condition of integrability 

distribution looks like a pure fluctuation distribution, with 

near zero requires b/a be positive. This leads us to iden- 

ill-defined mean value (b) the scaling for ||_B|| ^ I (mag- 

tify a necessary and sufficient condition for existence of a 

netic energy) is the same as that of 6 : P(||i?||) = \\B\p 

stationary dynamo as 

with J — b/ D — 1, where D is a "diffusion coefficient" for 
a > and ^ > DYNAMO (9) the finite-time Lyapunov exponent. On the contrary, the 
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Fig. 1. Result of the surrogate ID model: dtx — {b + (,{t))x — 
71^ with (C(t)C(i')> = 2a(5(t - t'). On the left side we show 
time series for a — 0.2, 7=1 and 3 different values of the 
parameter b. On the right side, the corresponding PDF and 
the theoretical curve corresponding to equation ([8}. 

bifurcation (II) is quite different in nature because of a 
well defined mean value for the magnetic field and fluc- 
tuations around this mean. Note also that if we consider 
the dynamo instability in absence of noise, a — 0, the 
two bifurcation threshold collapse. As stated by [5], the 
bifurcation corresponding to 5 > may be difficult to ob- 
serve in real experiments because, under the threshold, 
the presence of the Earth external magnetic field always 
gives rise to magnetic fluctuations qualitatively similar to 
that above the threshold (magnetic bursts separated by 
quiescent period). This effect can be taken into account 
by adding an additive noise to equation P]). 

In the theory of dynamical systems stability, the insta- 
bility criterion is usually associated with the existence of a 
positive Lyapunov exponent for the growth of the system 
energy: limlnB^ /2t = limlnS/t. It is possible to find 



this exponent, by multiplying equation ([5]) by Ini? and 
integrating with respect to B. This yields dt{\nB) — b, 
meaning that the Lyapunov exponent in our system is 
equal to b. The two instability criteria (existence of a nor- 
malizable solution or a positive Lyapunov exponent) are 
therefore identical provided a > 0, a necessary condition 
for integrability of the PDF at infinity. 

It is now interesting to discuss qualitatively the mean- 
ing of our main result ([9]). It is possible to show that for 
isotropic or axisymmetric velocity fiuctuations, the coeffi- 
cient a is positive. So we suspect that the main condition 
for existence of a dynamo is positivity of b. Therefore, in 
the limit of zero noise, the term proportional to fi is negli- 
gible and the dynamo threshold is only determined by the 
condition {dkVieieu) rj, > 0. Since the magnetic field will 
mainly grow in the direction given by the largest eigen- 
value of Sij — djVi and that the molecular diffusivity 
will tend to orientate the magnetic field along contract- 
ing direction, the dynamo threshold will be the same as 
in the deterministic case. Consider now a situation where 
you increase the noise level. Two different influences on 
the sign of b then result: one through the factor propor- 
tional by /i. According to the sign of this factor, it can 
therefore favor or hinder the dynamo. Another less obvi- 
ous influence is through vector orientation. Indeed, noise 
changes the distribution of magnetic field orientation. For 
example, if noise induces a flat distribution for e^, then 
{dkVieiek) rj, = Sii 0. Even more dramatic results can 
be obtained if the noise tends to align the vector along 
a direction of negative eigenvalue for Sij, since in that 
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case the factor {dkVieiek) ij, becomes negative, decreasing 
the dynamo threshold. From these two examples, one sees 
that in the presence of noise, the dynamo threshold can 
be completely disconnected from the dynamo threshold 
in the deterministic system ! Some caution is therefore in 
order when designing an experiment based on mean field 
measurements. 

Our approach gives a quantitative criterion on the dy- 
namo threshold (namely b > Q). However, its practical 
implementation requires the measure on the angular and 
position variables G(x, n, t) and the average of the tensor 
/i with this measure. One can obtain the measure by in- 
tegrating equation ([3]) with respect to B. Unfortunately 
the equation for G can not be solved in the general case 
and particular types of turbulence statistics have to be 
considered (isotropic, axisymmetric, etc.). Work is under 
progress to determine the angular measure in these simple 
cases. It is however interesting to note that this measure 
explicitly involves the tensors a and (3 defined in (U) and 
appearing in the Mean Field Equations ([5]). In that sense, 
the dynamo threshold depends on these tensors, albeit in a 
less explicit way than in the Mean Field Equations (MFE). 
It would therefore be interesting to confront threshold de- 
rived from (MFE), which are independent, and from our 
theory, to see what kind of error in the threshold determi- 
nation one can expect by using MFE instead of the true, 
non-perturbative theory. 
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